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Abstract 

Theoretical investigation of thermodynamic properties of electrically modulated monolayer 

graphene in the presence of a perpendicular magnetic field B is presented. The results obtained 
C3 ! 
'-*h ■ are compared with those of the conventional 2DEG. The one-dimentional periodic potential due 

to electric modulation lifts the degeneracy of the Landau Levels and converts them into bands 



ctf 



whose width oscillates as the function of B. We find Weiss type oscillations for small values of B 
and dHvA type oscillations at larger values values of B. These oscillations are more pronounced in 



Graphene than in conventional 2DEG system. These oscillations are less damped with temperature 
in Graphene compared with conventional 2DEG systems. 
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I. INTRODUCTION 

Graphene is a 2D-honeycomb lattice of carbon atoms. Its experimental realization has 
opened up new horizons in the field of condensed matter physics and material sciences. 
Unique electronic properties of Graphene make it substantially different from conventional 
2DEG systems. The quasi particles in graphene at low energies have a linear dispersion 
relation e^ = hvpk with characteristic velocity of vp = 10 6 m/s[l]. These quasi particles 
called massless Dirac Fermions, can be treated as electrons with zero mass or neutrinos 
having electronic charge. The zero mass property of charge carriers in graphene along 
with charge conjugation symmetry, results in many unusual transport phenomena such as 
anomalous Quantum Hall Effect, Chiral Tunneling and non-zero Berry's Phase[2-6]. The 
2D Dirac like spectrum was also confirmed recently by cyclotron resonance measurements 
in monolayer Graphene[l] and also by angle resolved photo electron spectroscopy [7]. 

Weiss oscillations [8, 9] appear in magnetoresistance when convential 2DEG is subjected 
to artificially created periodic potentials (either electric or magnetic) in submicron range. 
Electrical modulation can be carried out by depositing an array of parallel metallic strips on 
the surface[12] or through two interfering laser beams[13]. These Oscillations are the direct 
consequence of the commensurability of two different length scales namely the cyclotron orbit 



radius R c = y/2rrn e l 2 (where n e is the density of electrons, and / = ^JhjeB is the magnetic 
length) and the period of modulation a. Weiss oscillations occur in the small magnetic 
field range [10, 11] and are separate from dHvA(de Hass van Alphen) and SdH(Subnikov 
de Hass) type oscillations which occur at larger values of magnetic field. These oscillations 
are found to be periodic in the inverse magnetic field. It is interesting to study the effects 
of electrical modulation on Dirac electrons in graphene. Theoretical studies of tranport 
properties of Dirac electron in graphene subjected to electrical modulation were recently 
carried out and showed the appearance of Weiss oscillations in magnetoconductivity. In 
addition, the magnetoplasmon spectrum of monolayer graphene in the presence of electrical 
modulation was recently investigated [14]. 

In this work we investigate the effects of elecrical modulation on thermodynamic prop- 
erties of monolayer graphene and compare the results obtained with those of conventional 
2DEG system found in semiconductor hetrostructures. To this end, wehave determined the 
following thermodynamic quantities: The chemical potential, Helmholtz free energy, orbital 



magnetization , orbital magnetic susceptibility and electronic specific heat. The results are 
compared with those of the conventional 2DEG studied in[15] and[16]. 

This paper is arranged as follows. In section II, we give the formulation of the problem. 
The calculation of the thermodynamic quantities is given in section III and numerical results 
with discussion are presented in section IV. Finally the Conclusions are drawn at the end. 

II. FORMULATION 

We consider monolayer graphene in the xy —plane subjected to a magnetic field B along 
the z-direction. In the Landau guage, the unperturbed Dirac like Hamiltonian for single 
electron may be written as [6] 

H = v F a. (-ihV + eA) . (la) 

Here, a = {cr x ,a y } are the Pauli matrices and vp = 10 6 m/s characterizes the electron 
velocity, and A = (0, Bx, 0) is the vector potential. The normalized eigenfunctions of the 
Hamiltonian given in Eq.(l)[17? ] 

= e ikyy f-^ n _ x [( x + Xo )/l]\ 

nMy v / 2l7V 4> n [{x + x )/i] y { } 

where <\> n = expl -~ x L l H n (x), H n (x) are the Hermite Polynomials, L y is the normalization 

length in the y-direction, n is an integer corresponding to the Landau Level index and 
x = k y l 2 , is the center of the cyclotron orbit. The energy eigenvalues are 



v F HV2n r 
e n = = y/nfuj c (3) 

where tu c = ^f\i'^- is the cyclotron frequency of the Dirac electrons in graphene. To 
investigate the effects of modulation we write the Hamiltonian in the presence of modulation 
as 

H = H + U(x) (4) 

where U(x) is the one-dimensional periodic modulation potential along the x-axis and is 
given by 

U(x) = V cosKx. (5) 



K = — , a is the period of modulation and V is the constant modulation amplitude. To 
account for the weak modulation we take V to be an order of magnitude smaller than the 



Fermi Energy e° F = vpfikF, where &p = \/2im s is the magnitude of Fermi wave vector. Hence 
we can apply standard first order perturbation theory to determine the energy eigenvalues 
in the presence of modulation. The first order energy correction is 

£ n ,x = e n + U n cos Kx (6) 

Here, U n = ^exp(-|)[L n (x) + L n _i(x)], x = ^r and ' L n{x) and L n -i(x) are Laguerre 
polynomials. 

Although similar features in the energy spectrum have also been found in the 2DEG 
system under similar conditions, [16? J there are substantial differences between the two 
systems. Landau level spectrum of Dirac electrons depends on the square root of both 
magnetic field B and the Landau band index n against linear dependence in the case of 
conventional electronsin 2DEG. The energy eigenvalues in the presence of modulation given 
by Eq. (6) contains a term which is a linear combination of two succesive Laguerre polyno- 
mials with indices n and n — 1 , while conventional electrons obey a relation containg one 
Laguerre polynomial with index n. 

The modulation potential lifts the degeneracy of the Landau levels and broadens the for- 
merly sharpe levels into electric Landau bands. The electric modulation induced broadening 
of the energy spectrum is nonuniform. The Landau band width U n oscillates as a function 
of n since L n (x) is anoscillatory function of the index n. These landau bands become flat for 
different values of B. Flat bands occure for those values of B for which modulation strength 
becomes zero. By putting Un = one can get the flat band condition. 

exp(-|)[L n ( X ) + L n _ 1 (x)] = (7) 

using the asymptotic expression [17] 

ex P(- o ) L «(x) - i ■ cos(2 v ^ - -) (8) 

and L n (x) = L n _i(x), one obtains from Eqs (6) and (7) the following condition 

2i2 c = a(i-l/4), i = l,2,3, (9) 

where, R c = kpl 2 , is the classical cyclotron orbit. From Eqs (6) and (8) it can be observed 
that, in the large n limit electron bandwidth oscillates sinosoidally and is periodic in 1/B, 



for fixed values of n and a. When n is small bandwidth still oscillates, but the condition (9) 
no longer holds because neigther eq. (8) nor L n (x) — L n _i(x) is valid. Interestingly, for low 
values of B, when many Landau levels are filled, both the systems have the same flat band 
condition [15]. 

It is well known that in the absence of modulation the density of states (DOS) consists 
of a series of delta functions at energies equal to e n . The addition of a weak periodic electric 
modulation however modifies the former delta functions leading to DOS broadening . The 
density of states D[e) are given by [18] 
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where, 9{x) is a unit Heaviside step function and A is the area of the sample. 



III. EQUILIBRIUM THERMODYNAMIC QUANTITIES 



We have determined the electronic contribution to the equillibrium thermodynamic prop- 
erties of monolayer graphene subjected to a perpendicular magnetic field and weak electric 
modulation. The thermodynamic quantities calculated are chemical potential, Helmholtz 
free energy, electronic specific heat, orbital magnetization and orbital magnetic susceptibil- 
ity. 

The magnetid field (B) and temperature (T) dependent chemical potential \x = fi(B,T) 
of a system can be determined by inverting the following relation 



oo 
N= JD(e)f{e)de 



(11) 



where the Fermi Dirac distribution function / (e) is 

' e — \i 



/(e) 



cxp 



k B T 



(12) 



ks is the Boltzmann's constant and N is the total number of electrons. Hence change in 
the D(e) will affect fi(B,T). Substituting Eq.(9) into Eq.(ll) we obtain 



i 



dx 






1 + x n exp[z n x]) 1 



(13) 



Here Xn = exp 



k B T 



and z n = \U n \ /(kBT). Eq.(12) can be used for both modulated and 
unmodulated systems (z n = 0). We solve this equation numerically in order to obtain the 
chemical potential fi(B,T). We are able to determine the change in the chemical potential 
due to the electric modulation. Once the chemical potential and the density of states are 
known, the free energy F of the system can be calculated. From there on the thermodynamic 
properties of the system can be obtained from the free energy by taking the appropriate 
derivatives. For a system of non-interacting fermions, the Helmholtz free energy is given by 
[19] 



F = fiN-k B T /*£>(£) In l + exp(^-\ 

o 



de (14) 



The density of states D(e) is the central quantity in the above expression. The expression 
for D(e) in graphene is different from that in conventional 2DEG due to the difference in the 
energy spectrum in the two cases. This difference will affect the electronic contribution in 
the thermodynamic properties in the two systems determined from the following free energy 
for the electrically modulated graphene system 

i 
F = fiN- k B T —Y, / f, f In [1 + Xn 1 exp (~z n x)] (15) 



n=0" l 



From Eq. 15, one can calculate the electronic comtribution of magnitization for both 
graphen and 2DEG systems as M = — (§§) A N and specific heat as C v = —T ( |^§ j 
The electronic contribution to susceptibility is obtained directly from x — ~ (d 2 F/dB 2 ) . 

IV. RESULTS AND DISCUSSION 

Numerical study of thermodynamic properties for monolayer graphene system subjected 
to electrical modulation is presented. We have also plotted the same quantities for the 
2DEG system This is to facilitate comparison and was also a check on our numerical 
program. For the 2DEG parameters for GaAs are used. We have taken n s = 3.16 x 10 15 m -2 
and a = 382nm. For electrical modulation we have taken Vo = lmeV. Thus our 2DEG 
results are those already given in[16? ]. Modulation induced effects on thermodynamic 
quantities can be highlighted by calculating the difference between the modulated case and 
the unmodulated case in each system. 
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In Figures 1 — 5 we have plotted the change in various thermodynamic properties due to 
electric potential at temperatures of T = 2K (full curve) and T = 6K (broken curve). for both 
conventional 2DEG system and graphene system. These figures were scaled to approperiate 
values to make them appear dimensionless. 

In Fig.(l), we have plotted the change in chemical potential versus magnetic field at 
temperatures 2fT(straight) and 6i^(broken). For Conventional 2DEG system for B < 0.3T 
oscillations depend very weakly on temperature, which is a clear signature of Weiss type 
Oscillations. Where as for B > 0.3T, the oscillations depends strongly on temperature, in 
particular they die out at 6K, a clear signature of dHvA type oscillations. Furthermore, 
the zeros in the chemical potential are in close agreement as predicted by the flat band 
condition Eq.(9). A similar behavior is expected for Graphene system. But for Graphene 
the value of B defining the boundary between the two oscillatory phenomena is quite low 
(It lies some where between O.land 0.15T). For smaller values of B Weiss type oscillations 
are present and the amplitude of the oscillations remain essentially the same at different 
temperatures. For larger values of B, the familiar dHvA-type oscillations are present ,as 
the amplitude of oscillations is reduced considerably at comparatively higher temperature. 
However the oscillatory phenomenon still persists, contrary to the conventional 2D system in 
which oscillations completely die out at 6K. In comparison we can say, Graphene system is 
more sensitive to the magnetic field and less sensitive to temperature, than the conventional 
2DEG system. This difference arises mainly due to the difference in the Landau level energies 
of the two systems and due to the presence of an additional Laguerre Polynomial term in 
the modulation contribution to the energy spectrum for Graphene system 

The Free energy is shown in Fig. (2), for the two systems. To make y-axis dimensionless, 
Free energy has been scaled using Fq = ^NEp. It can be seen that at small values of B 
periodically modulated potential induces temperature independent Weiss type oscillations, 
with zeros occurring at their respective flat band conditions. Weiss Oscillations are more 
prounounced in Graphene system, significantly the amplitude of Weiss oscillations for the 
graphene system remains unchanged at higher temperature, contrary to the 2DEG in which 
damping may be observed. The familiar dHvA type oscillations are observed for higher 
values of B. As in the case of the Chemical Potential, again the dHvA type oscillations 
starts quite early. The first period for the dHvA type oscillations starts from B = 0.3T and 
extends up to 0.6T for the standard 2DEG system, while for graphene the first period of 



dHvA type starts near B = 0.175T and terminates at 0.27T. 

In Figs. (3) and (4) we have plotted the changes in the magnetization AM and the suscep- 
tibility A% against the magnetic field. Both the quantities has been approperiately scaled 
to appear dimensionless. At low B oscillations having their origin in the commensurability 
of two length scales, and are effected weakly by temperature, having zeros as given by their 
respective flat band conditions. For higher values of B, dHvA oscillations are present at 
lower temprature (2K), with amplitude becomeing zero for the conventional 2DEG system 
while reduced considerably for the Graphene system at higher temprature (QK). 

In fig. 5 we plot change in the electronic specific heat capacity against magnetic field. 
y-axis has been scaled using C e \ = Nks, to appear dimensionless. In both systems the 
Weiss type oscillations are not large effects, however the damping behavior of dHvA type 
oscillations is clearly observeable. 

V. CONCLUSIONS 

We have presented a study of the thermodynamic properties of monolayer graphene sys- 
tem and compared the results with those of the conventional 2DEG. The commensurability 
oscillations(Weiss type) and dHvA type oscillations are reflected in all the thermodynamic 
quantities under consideration in this work for the two systems. However, these effects are 
more prounounced in graphene system in the sense that the oscillations in the thermody- 
namic quantities are more robust against temperature. We can therefore say that Graphene 
system is less sensitive to temprature and more sensitive to the magnetic field. This differnce 
arises because of the different nature of the quasiparticles in the two systems. 

VI. REFERENCES 



[1] R. S. Deacon, K. C Chuang, R. J. Nicholas K. S. Novoselov A. K. Geim. Phy. Rev. B 76, 

081406 (R) (2007). 
[2] Novoselov K S, Geim A K, Morozov S V , Jiang D, Katsnelson M I, Grigorieva I V, Dubonson 

S V and Firsov A A 2005 Nature 438,197 - 200. 

8 



[3] Y. Zhang, Y. W. Tan, H.L. Stormer, and P. Kim, ibid. 438, 201 (2005). 

[4] Y. Zhang and T. Ando, Phys. Rev.. B 65, 245420 (2002). 

[5] V. P. Guysynin and S. G. Sharapov, Phy. Rev. Lett. 95, 146801 (2005). 

[6] M. I. Kattsenelson, K. S. Novoselov and A. K. Geim Nature Physics 2, 620 - 625 (2006). 

[7] S. Y. Zhou, G.-H. Gweon, J. Graf, A. V. Fedrov, C. D. Spataru, R. D. Diehl, Y. Kopolevich, 

D. H. Lee, S. G. Louie, and A. Lanzara, Nat. Mater. 2, 595 (2006). 
[8] Gerhardts R R, Weiss D and von Klitzing K 1989 Phys. Rev. Lett. 62 1173 
[9] Carmona H A, Geim A K, Nogaret A, Main P C, Foster T J, Henini M, Beaumont S P and 
Blamire M G 1995 Phys. Rev. Lett. 74 3009 
[10] Weiss D, von Klitzing K, Ploog K and Weimann G 1989 Europhys. Lett. 8 179. 
[11] Weiss D, von Klitzing K, Ploog K and Weimann G 1989 High Magnetic Fields in Semiconduc- 
tor Physics II (Springer Series in Solid-State Sciences 87) ed G Landwehr (Berlin: Springer) 
p 357 
[12] D. Weiss, K. von Klitzing, K. Ploog, and G. Weinmmann, Euro-phys. Lett. 8, 179 (1989) 
[13] R. W. Winkler, J. P. Kotthaus, and K. Ploog, Phy. Rev. Lett. 62, 1177 (1989) 
[14] A. Mataulis and F. M. Peeters, Phy. Rev. B 75,1254929 (2007) 
[15] S M Stewart and Chao Zhang J. Phy: Condens. Matter 10 (1998) 5545-5566 
[16] F. M. Peeters and P. Vasilopoulous, Phy. Rev. B 46, 4667 (1992) 

[17] Tables of Integrals, Series, and Products, edited by I. S. Gradshteyn and I. M. Ryzhik (Aca- 
demic, San Diego, 1980) 
[18] M Tahir, K Sabeeh and A Mackinnon, J. Phys.: Condens. Matter 19, 406226(2007) 
[19] Patheria R K 1972 Statistical Mechanics (Oxford: Pergamon) p 215 



